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ABSTRACT 

We study the virialization of the cosmic structures in the framework of fiat 
cosmological models where the dark energy component plays an important role 
in the global dynamics of the universe. In particular, our analysis focuses on 
the study of the spherical matter perturbations, as the latter decouple from 
the background expansion and start to "turn around" and finally collapse. 
We generalize this procedure taking into account models with an equation 
of state which vary with time, and provide a complete formulation of the 
cluster virialization attempting to address the nonlinear regime of structure 
formation. In particular, assuming that clusters have collapsed prior to the 
epoch of Zf ~ 1.4, in which the most distant cluster has been found, we show 
that the behavior of the spherical collapse model depends on the functional 
form of the equation of state. 
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1 INTRODUCTION 

It is well known that the available high quality cosmo- 
logical data (Type la supernovae, CMB, etc.) are well 
fitted by an emerging "standard model" , which contains 
cold dark matter (CDM) to explain clustering and an 
extra component with negative pressure, the dark en- 
ergy, to explain the observed accelerated cosmic expan- 
sion. The last decades there have been many theoretical 
speculations regarding the nature of the above exotic 
dark energy. Most of the authors claim that a scalar field 
which rolls down the potential V{4>) (Ratra & Peebles 
1988; Caldwell, Dave, Steinhardt 1998; Peebles & Ra- 
tra 2003) could resembles that of the dark energy. The 
cosmic equation of state is given by pq = wpq, where 
pg is the pressure and pq is the corresponding density 
of the dark energy. Owing to the fact that w < — 1/3 is 
related to the potential of the dark energy field, having 
an indication about its value may help us to understand 
the nature of the dark energy. 

However, a serious issue here is how the large scale 
structures and in particular galaxy clusters form. The 
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cluster distribution traces basically scales that have not 
yet undergone the non-linear phase of gravitationally 
clustering, thus simplifying their connections to the ini- 
tial conditions of cosmic structure formation. The so 
called spherical collapse model, which has a long history 
in cosmology, is a simple but still a fundamental tool to 
describe the growth of bound systems in the universe 
via gravitation instability (Gunn & Gott 1972). In the 
last decade many authors have been involved in this 
kind of studies and have found that the main param- 
eters of the spherical collapse model such as the ratio 
between the final and the turn around radius (hereafter 
collapse factor), is affected by the dark energy (Lahav 
et al. 1991; Wang & Steinhardt 1998; Iliev & Shapiro 
2001; Battye & Weller 2003; Basilakos 2003; Weinberg & 
Kamionkowski 2003; Mota & van de Bruck 2004; Horel- 
lou & Berge 2005; Zeng & Gao 2005; Maor & Lahav 
2005; Percival 2005; Wang 2006; Nunes & Mota 2006). 

The aim of this work is along the same lines, at- 
tempting to investigate the cluster formation processes 
by generalizing the non-linear spherical model for a fam- 
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ily of cosmological models with an equation of state 
parameter being a function of time, w — w(t). This 
can help us to understand better the theoretical ex- 
pectations of negative pressure models as well as the 
variants from the Quintessence (w = wo = const with 
— 1 < wo < — 1/3) and Phantom (w = wo = const with 
wq < — 1) case respectively. 

The structure of the paper is as follows. The basic 
elements of the cosmological equations are presented in 
section 2. Sections 3 and 4 outline the spherical collapse 
analysis in models where w is a function of time and 
finally, we draw our conclusions in section 5. 



2 THE BASIC COSMOLOGICAL 
EQUATIONS 

For homogeneous and isotropic cosmologies, driven by 
non relativistic matter and an exotic fluid with equation 
of state, pq — w(a)pQ with pq < 0, the Einstein field 
equations can be given by: 

2 8ttG, n k 



(1) 



with k = —1,0 or 1 for open, flat and closed universe 
respectively and 



a 
a 



-4ttG 



(«(a) + I) 



PQ + gPn 



(2) 



where a(t) is the scale factor, p m = pmooT 3, is the back- 
ground matter density and pq = PQof(ct) is the dark 
energy density, with: 



f(a) = exp 



1 + w(u) 



du 



(3) 



Thus, the scale factor evolves according to Fried- 
mann equation: H 2 = (a/a) 2 . The relation between the 
time and the scale factor is given by 

* = - J— (4) 

da H(a)a ' V ' 

where the Hubble parameter is written: H(a) = 
HoE(a) and Ho is the Hubble constant with 

E{a) = [n m a- 3 + Q k a- 2 + n Q f(a)] 1/2 , (5) 

while Qm = 8-rrGpmo/SHo (density parameter), Q, k = 
—k/Ho (curvature parameter), Q,q = 8ttGpqo/3Ho 
(dark energy parameter) at the present time with Q m + 
fifc + Qq = 1. 

In addition, to Q m (a) also Qq (a) could evolve with 
the scale factor as 



n m (a) = 



and flQ(a) = 



(6) 



E 2 (a) wy ' E 2 (a) ' 

Note that in this paper we consider a spatially flat 
(k — 0) low-f2 m = 0.3 cosmology with Ho = 72 Km 



s _1 Mpc _1 which is in agreement with the cosmologi- 
cal parameters found from the recent observations (see 
Freedman et al. 2001; Spergel et al. 2003; Tegmark et 
al. 2004; Basilakos & Plionis 2005; Spergel et al. 2006 
and references therein). Finally, in order to address the 
negative pressure term it is essential to define the func- 
tional form of the equation parameter w — w(a) [see 
section 4]. 



3 VIRILIZATION IN THE SPHERICAL 
MODEL 

The spherical collapse model is still a powerful tool, de- 
spite its simplicity, for understanding how a small spher- 
ical patch of homogeneous overdensity forms a bound 
system via gravitation instability. Technically speaking, 
the basic cosmological equations, mentioned before, are 
valid either for the entire universe or for homogeneous 
spherical perturbations [by replacing the scale factor 
with radius R(t)] 



-4nG 



1\ 1 

W ( R ) + q ) PQc + oPmc 



(7) 



R 

R LV~'"" y ' 37 ''^ ' 3 r 

where p mc and pq c is the time-varying matter and 
dark energy density respectively (for spherical pertur- 
bations) . 

We study the cluster virialization in models with 
dark energy, generalizing the notations of Lahav et al. 
(1991), Wang & Steinhardt (1998), Mainini et al. (2003); 
Basilakos (2003), Lokas (2003) Mota & van de Bruck 

(2004) , Horellou & Berge (2005), Zeng & Gao (2005), 
Maor & Lahav (2005), Bartelmann, Doran & Wetterich 

(2005) , Nunes & Mota (2006) and Wang (2006) in order 
to take into account models with a time varying equa- 
tion of state. Thus, in this section, we review only some 
basic concepts of the problem. Within the framework 
of the spherical collapse model we assume a spherical 
mass overdensity shell, utilizing both the virial theo- 
rem T c — —\Ugc + Uq c and the energy conservation 
T c j + Uocj + UqcJ = Ucc,t + UQc,t where, T c is the ki- 
netic energy, Ugc = —3GM 2 /5R is the potential energy 
and Uq c is the potential energy associated with the dark 
energy for the spherical overdensity. In particular, the 
potential energy induced by the dark energy component 
(see Horellou & Berge 2005 and references therein) is 



Uqc = (1 + 3w)p Qc — — — R 



(8) 



Using the above formulation we can obtain a cubic equa- 
tion which relates the ratio between the final (virial) Rf 
and the turn-around outer radius R t : 



2mX 3 - (2 + n 2 )A + 1 = 



(9) 
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where A = ^ is the collapse factor, 



ni — — (3w + 1) 



and 



n 2 



with 



pmc.t 
pm,t 



£fi m a t 

-Rt N ~ 3 
at . 



(10) 



(11) 



(12) 



(the definition of the a t and at factors are presented 
bellow). 

However, we would like to point out that there is 
some confusion in the literature regarding eq.(9) which 
is based on energy conservation. Indeed recently it has 
been shown (Maor & Lahav 2005) that the above for- 
mulation is problematic due to the fact that the total 
energy of the bound system is not conserved in dark 
energy models with w 7^ — 1. In order to avoid these 
systematic effects in this work we utilize two different as- 
sumptions. First of all, we assume that when the matter 
epoch just begins, the dark energy moves synchronously 
with ordinary matter on both, the Hubble scale and the 
galaxy cluster scale, the so called clustered dark energy 
scenario (see Zeng & Gao 2005 and references therein). 
It is interesting to mention that the synchronous dark 
energy evolution (hereafter clustered) was designed to 
conserve energy. In particular, Maor & Lahav (2005) ad- 
dress the issue of the clustered dark energy model based 
on the following assumptions: (i) clustered quintessence 
considering that the whole system virializes (matter and 
dark energy) and (ii) the dark energy remains clustered 
but now only the matter virializes (for more details see 
their section 4). Note, that in this work we are utilizing 
the second possibility. In that case the equation which 
defines the collapse factor becomes: 



(13) 



(14) 



(l + g)A-|(l-3HA- 3tu = i 

where the factor q is given by 
q _ PQ.t _ PQ,t _ v_ 

pmc.t CPm.t C 

(the definition of the v parameter is presented in section 
3.1). 

If we assume that the dark energy component in 
a galaxy cluster scale can be treated as being homoge- 
neous then we get the following cubic equation defined 
by Maor & Lahav (2005) 



2qA 3 -(l + g)A+-=0 , 



(15) 



Although it was pointed out recently (see Wang 2006) 
that if the value of q is less than 0.01 then the prob- 
lem of energy conservation does not really affect the 



virialization process and thus eq.(9) is still a good ap- 
proximation. 

From now on, we will call at the scale factor of the 
universe where the overdensity reaches at its maximum 
expansion and at the scale factor in which the sphere 
virializes, while R t and R{ the corresponding radii of 
the spherical overdensity. Note that p mc ,t is the matter 
density which contains the spherical overdensity at the 
turn around time while p mjt is the background matter 
density at the same epoch. Therefore, for bound pertur- 
bations which do not expand forever, the time needed 
to re-collapse is twice the turn-around time, tt = 2t t . 
Doing so and taking into account eq.(4), it is easy to 
estimate the relation between the at and a% 



1 



-da = 2 



1 



-da 



H(a)a- "J Q H(a)a~~ ' (16) 

In order to reduce the parameter space of the over- 
all problem we put further observational constrains. 
Indeed, at can be defined from the literature assum- 
ing that clusters have collapsed prior to the epoch of 
Zf ~ 1.4 in which the most distant cluster has been 
found (Mullis et al. 2005; Stanford et al. 2006). There- 
fore, considering at ~ 0.41 and utilizing eq.(16), it is 
routine to obtain the scale factor at the turn around 
time*. In the case of a A cosmology we get an analyti- 
cal solution: 

sin~ 1 h (^a A / 2 ^/vo^j = 2sin~ 1 h ^a 3 ^ 2 ^oj (17) 

where vq = (1 — f2 m )/fi m . While for the general problem 
we have to solve equation (16) numerically. The ratio 
between the scale factors converges to the Einstein de 
Sitter value = 2 2/3 at high redshifts. 



3.1 The rescaled equations 

In this section we present some of the basic concepts of 
the spherical collapse model. In particular, we assume 
a spherical overdensity which contains a dark energy 
component which behaves either as clustered or homo- 
geneous. For a flat cosmology (k = 0), using the basic 
differential equations (see eq.l and 7) and performing 
the following transformations 



a , R 

x = — and y = — 

at Rt 



(18) 



the evolution of the scale factors of the background and 
of the perturbation are governed respectively by the fol- 
lowing two equations: 



x 2 = HtQ, m ,t [Q. m {x)x] 1 



The epoch of the turn around is roughly zt ~ 2.5 



(19) 
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and 



with 



y = - 



H\ Qm,t 



+ vyi{x,y) 



(20) 



f R Y- 

pmc — Pmc.t yj£ J — 



Om,t 



and 

I{x,y) 

with 



[l + 3w(R(y))]^l 
[l + 3w(x)]f(x) 



PQ,t 



1 + vx z f(x) 



y 2 = -tf t 2 fi m , t [-(y- 1 + G(y) - C] 



where C is the integration constant and 



G(y) 



/( 



«t)/ 



y[l + 3w(R(y))]f(R(y))dy 



(2)' 



y 2 /(^(y)) 

/(«t) 



(28) 



Note, that in order to obtain the above set of equa- 
tions we have used the following relations: 



(21) 

Clustered DE 
Homogeneous DE 



(22) 



Finally, below we present the general integral equa- 
tion which governs the behavior of the density contrast 
C at the turn around epoch, for generic flat cosmological 
models, as a function of the equation of state (which also 
depends on time) and the perturbation collapse time 

1/2 

An, 

(29) 



f 1 

Jo 



(23) 

Pm,t S'm,t 

The function R(y) is given by the combination of cq.(12) 
and eq.(18): R(y) = R t y = C 1/3 a t JA Finally, fJ m (x) is 
given by 

1 



(24) 



3.2 The general solution for the clustered 
dark energy scenario 

In this section we present our analytical solution of the 
£ parameter by integrating the above system of differen- 
tial equations, (eq.19 and eq.20), using at the same time 
the boundary conditions of (dy/dx) x= i = and y = 1. 
In the case of a homogeneous dark energy the above 
system is solved only numerically. The novelty of our 
approach here is that for the clustered dark energy case 
the system can be solved analytically including models, 
where the dark energy parameter w is a function of the 
cosmic time. Indeed, due to the fact that now the second 
part of eq.(20) is a function only of y we can perform 
easily the integration 



(25) 



Considering now that the functional form of f{R(y)) is 
exponential (see the appendix) we have the following 
useful formula: 

/ y[l + 3w(R(y))]f(R(y))dy = -2 / yf(R(y))dy 

-fy 2d -^dy = -y 2 W(y)) ■ ( 3 

Therefore, the basic differential equation for the evolu- 
tion of the overdensity perturbations under the frame- 
work of the boundary conditions [described before C = 
£ + P{l)v\ takes the form 

'Ay\ 2 (y- 1 + vP(y) - ( - P(l)v 



[xdm(x)] 1 



(27) 



= Jl [xtt m (x) 



dy 

x ' 2 dx 



4 SPECIFIC DARK ENERGY MODELS 

Knowing the functional form of the function f(x), the 
collapse scale factor (in our case at ~ 0.41), the turn 
around scale factor a t (from eq. 16) and assuming a 
low matter density flat cosmological model with 57 m = 
1 — Oq = 0.3, we can obtain the parameter £ solving 
numerically the equation (29). Note that in this work, 
we deal with four different type of equations of state, 
which we present in the coming sections. 



4.1 Quintessence - Phantom models 

In this case the equation of state is constant (see for a 
review, Peebles & Ratra 2003; Caldwell 2002; Corasan- 
iti et al. 2004). If — 1 < wo < — 1/3 we have the so called 
Quintessence models while for wo < —1 we get the so 
called Phantom cosmologies (hereafter QP-models). De- 
spite the fact that for these models the cluster virializa- 
tion has been investigated thoroughly in several papers 
we have decided to re-estimate it in order to understand 
the variants from the w — w(t) case. In particular, for 
w(a) — wo we get f(a) = a _3 ( 1 + IU0 ' and the overall 
problem (eq. 29) reduces to the integral equation found 
by Zeng & Gao (2005). In particular, the basic factors 
become 



P(y) = f+vo y -3w -i 
and 

i 

Qm(x) = 



(30) 



(31) 



1 + i>x- 3w ° 

In Figure 1 we present for the clustered (solid line) and 
homogeneous dark energy respectively (dashed line) the 
£ solution (top panel), as a function of the equation of 
state parameter and it is obvious that the two models 
give almost the same results^. Also Figure 1 (bottom 

t For wo < —0.9 the ( solution tends to the Einstein-de 
Sitter case C > (^) 2 
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Figure 1. The upper panel shows the f parameter as a func- 
tion of wo- Note that the solid line presents the clustered 
dark energy scenario while the dashed line corresponds to 
the homogeneous dark energy. The lower panel shows the 
corresponding density contrast at virialization. for a flat cos- 
mology, Qm = 1 — Qq = 0.3 with a constant equation of 
state. Finally, the insert plot shows the behavior of the q pa- 
rameter for different wo- Finally, the index P and Q reads 
Phantom and Quintessence respectively. 



panel) shows the dependence of the density contrast at 
virialization on w. 

p mc , c (32) 



A vir 



pm,f A 3 

where p mo ,f is the matter density in the virialized struc- 
ture while p m _f is the background matter density at 
the same epoch. The factor A v i r is a key parameter 
in this kind of studies because we can compare the 
predictions provided by the spherical collapse model 
with observations. In this framework, it is also inter- 
esting to note that the collapse factor is in between 
0.48 < A < 0.50 while for the homogeneous dark en- 
ergy we get 0.47 < A < 0.50. Both cases tend to the 
standard value 0.5 in an Einstein de Sitter universe in 
agreement with previous studies (Maor & Lahav 2005; 
Wang 2006). This is to be expected simply because at 
large redshifts matter dominates over the dark energy 
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Figure 2. The non linear infall velocity in the virialized 
regime as a function of the equation of state parameter. 



in the universe which means that the parameter q has 
small values (see the insert plot of Fig. 1). Note that 
A-models can be described by quintessence models with 
wo strictly equal to -1 and thus, eq.(29) is written: 

1/2 

dy 



V 



( + vy 3 -(( + v)y 



ln(Vl + ^ + -Jv) 



2/3 



(33) 



In order to explore further the virialized structures 
we investigate the connection between the infall velocity 
at the virialized epoch and the equation of state. It is 
well known from the literature that the non linear infall 
velocity field is described by the following expression 
(see Lilje & Efstathiou 1989; Croft, Dalton & Efstathiou 
1999): 



V 1 



5(r) 



(34) 



3""'"' uv "'' [1 + <5(r)] - 25 ' 

where S(r) is the fluctuation field. It should be noted 
that in this work we have used the generic expression for 
the growth index (~ f2^) defined by Wang & Steinhardt 
(1998) with 







- + 



3(l-wo)(l-3w /2) 



(l-n m ) .(35) 



5 - w /(l - wo) ' 125(1 - 6w /5) 3 
Due to interplay between the infall velocity and the 
virial radius, in Figure 2 we present the correspond- 
ing ratio between them at the virialized epoch, [r = Rf, 
a = Of and 8(R{) ~ A v j r ], as a function of the equation 
of state parameter. Therefore, it becomes evident that 
for models with wo < — 1 , the so called Phantom models 
(see Caldwell 2002), the infall velocity is unaffected re- 
ally by the equation of state. When we use models with 
w > — 1 the infall velocity is a decreasing function of w 
(the difference is ~ 37%). The latter is to be expected 
because of the A v i r behavior (see Figure 1). 
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4.2 Models with time variation w 

The last few years there have been many theoretical 
speculations regarding the nature of the exotic "dark 
energy" . Various candidates have been proposed in the 
literature, among which a dynamical scalar field act- 
ing as vacuum energy (Ozer & Taha 1987). Under this 
framework, high energy field theories generically show 
that the equation of state of such a dark energy is a 
function of the cosmic time. If this is the case then the 
basic equations regarding the cluster virilization in mod- 
els with a time varying equation of state become more 
complicated than in models with constant w. Note that 
in the literature due to the absence of a physically well- 
motivated fundamental theory, there is plenty of dark 
energy models that can fit current observational data 
(for a review see Liberato & Rosenfeld 2006). In this 
work as an example, we use several parameterizations 
regarding the dark energy component in order to evalu- 
ate the integrals of eq.(3) and eq. (29). In particular, we 
consider a simple expression for the equation of state: 

w(a) = wo + wiX(a) 

which is available in the literature (see for example Go- 
liath et al. 2001; Linder 2003; Cepa 2004; Liberato & 
Rosenfeld 2006 and references therein). The function 
X(a) satisfies the following constrain at the present 
epoch: X(l) — and thus, wo < —1/3. Here we sam- 
ple the parameters as follows: w G [—1.05, —0.35] and 
w-i G [0.1, 1.5] in steps of 0.1. To this end, we need to 
emphasize that our approach is powerful in a sense that 
for different parametric forms of w(a) the system of 
eq. (19) and eq. (20) can be solved analytically in the 
framework of a clustered dark energy. Finally, in this 
work wc assume that the parameters (wo , wi ) and the 
functional form of the equation of state parameter re- 
main the same either for the entire universe or for the 
spherical perturbations. 




0.5 1 



Figure 3. Possible pairs (wq, u>i) defining the collapse factor 
A for the clustered case, based on w(a) = uio + m>i(1 — a) [see 
section 4.2.1], The lines correspond to wo = —0.35 (solid), 
wo = —0.45 (short dashed), wo = —0.55 (dashed) and wo = 
-0.65 (dot dashed). 

For the homogeneous dark energy scenario we find 
that the collapse factor lies in the range 0.47 < A < 0.50 
while for the clustered case we get 0.30 < A < 0.50. 
In particular, Fig. 3 shows the behavior of the col- 
lapse factor for possible pairs of (w ,wi) starting from 
wo = —0.35 (solid line), wo = —0.45 (short dashed 
line), wo = —0.55 (dashed line) and wo = —0.65 (dot 
dashed line). From the figure, it becomes evident that 
for large values of wi the collapse factor is significant less 
than 0.5 which means that under the framework of the 
present equation of state, a candidate structure (large 
scale overdensity) is located in a large density regime 
(A v i r > 400) and at the end the tendency is to collapse 
in a more bound system, with respect to the QP-models 
(w = const). Additionally, it is interesting to mention 
that for small values of wi < 0.55 the structures re- 
veal a similar dynamical behavior, 0.48 < A < 0.50, as 
described in section 4.1. 



4-2.1 The w(a) — wo + wi(l — a) model 



4.2.2 The w(a) = w + wi(a 1 — 1) model 



In this case, (see Linder 2003; Cepa 2004) where wo and 
wi are constants, using eq.(3), eq.(24) and eq.(28) we 
derive the following basic functions of the differential 
equation eq.(29): 

^(q.) _ a -3(l+u> +u>i) e 3«;i(a-l) 
1 



fi m (a;) 



and 



P{y) 



(36) 
(37) 



c 1+w 



V 



3(1110+™!)+! 



(38) 



Using now the equation of state derived by Goliath et 
al. (2001) then the formulas described before become: 



jV a ) _ Q ,-3(l+u>o-u>l) e 3iDi(c« 



and 



P(y) 



1 



1 + ifx~ 3 ( w o- w i)e 3wi ( x 



£l+MJ0-wig3uJia t 1 (C 1 ' 3 !/ * — 1) 



y 



3(w — u>i) + l 



(39) 
(40) 

(41) 



Following the same paradigm as in the previous case in 
Fig.4 we present for the clustered case the distribution 
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Figure 4. The collapse factor A, based on w(a) = wq + 
wi(a^ 1 — 1) [see section 4.2.2]. The lines correspond to the 
same values of wq as in figure 3. 



of A. It is obvious that also this equation of state affects 
the virialization of the large scale structures in the same 
way as before (see section 4.2.1), but the corresponding 
wi parameter takes much lower values. For wi < 0.15 
the collapse factor tends to 0.5. Note, that for the homo- 
geneous dark energy the collapse factor is in the range 
0.45 < A < 0.49. 



rise to the observed large-scale structure of the universe. 
In this work assuming that the dark energy moves syn- 
chronously with ordinary matter on both the Hubble 
scale and the galaxy cluster scale (clustered dark en- 
ergy) , we have treated analytically the nonlinear spher- 
ical collapse scenario considering different models with 
negative pressure which contains also a time varying 
equation of state. We verify that having an overdense 
region of matter, which will create a cluster of galax- 
ies prior to the epoch of Zf ~ 1.4, dark energy affects 
the virialization of the large scale structures in the fol- 
lowing manner: flat cosmological models that contain 
either (i) a constant equation state or (ii) an equation 
of state with w = w + ^Wi(l — a)a the virialization ra- 
dius divided by the turn around radius is in between 
0.47 < R{/Rt < 0.50 and the density contrast at virial- 
ization lies in the interval 171 < A v i r < 240. It is inter- 
esting to mention that the same behavior is also found 
for the homogeneous dark energy. Finally, we find that 
the following equations of state: w = wo + wi(l — a) 
and w = wo + wi(a _1 — 1) correspond to relative 
larger A v i r 's and thus, produce more bound systems 
(0.30 < Rf/Rt < 0.50) with respect to the other two 
equations of state mentioned before. While for the ho- 
mogeneous dark energy we found structures where the 
collapse factor is in between 0.45 < Rf/Rt < 0.49. 



4-2.3 The w(a) — wo + wi(l — a)a model 

In this case, (see Liberato & Rosenfeld 2006) we get: 

(42) 



f(a)=a- S( - 1+ ^e^^ 



O m (x) = 



and 



1 + ra- 3t »«eT 1 ( I - 1 » 2 



(43) 



= Q ~ ^ • (44) 

Performing once again the same constrains as before 
the virialization radius to the turn around radius lies in 
the range 0.46 < A < 0.50. This result is in agreement 
with those derived in the QP-models, despite the fact 
that the two dark energy models have completely dif- 
ferent equations of state. Finally, considering now the 
homogeneous dark energy scenario the collapse factor is 
almost in the same interval i.e, 0.45 < A < 0.48. 



5 CONCLUSIONS 

The launch of the recent observational data has brought 
great progress in understanding the physics of gravita- 
tional collapse as well as the mechanisms that have given 
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APPENDIX A 

Without wanting to appear too pedagogical, we remind 
the reader of some basic elements of Algebra. Given a 
cubic equation: A 3 + aiA 2 + a- 2 \ + a-j, — 0, let D be the 
discriminant: 



T> — a\a\ — 4a?> — Aa\ a 3 — 27ci3 + 18ai<Z2a3 



(45) 



and 



3 9 27 
Xi = —a 1 + -aia 2 —a 3 , 



x 2 = - - 



3V3t> 



If V > 0, all roots are real (irreducible case). In that 
case Ai , A2 and A3 can be written: 
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01 2r 1 / 3 



9- {fi- 1)tt 



M = 1,2, 3(46) 



where r = y^cf + an d # — cos _1 (xi/r). 

In this study, we derive analytically the exact solu- 
tion of the basic cubic equation (eq.15), having poly- 
nomial parameters: oi = 0, ai — —(1 + q)/2q and 
a-i = l/4<?. Then the discriminant of eq.(15) is: 

8(l + q) 3 -27g 2 



V(q) 



16g 3 



Of course, in order to obtain physically acceptance re- 
sults we need to take q > which gives T>(q) > 0. 
Therefore, all roots of the cubic equation are real (irre- 
ducible case) but one of them A3 corresponds to expand- 
ing shells. It is obvious that for q — > the above solu- 
tion tends to the Einstein-de Sitter case (A3 — > 0.50), 
as it should. 



APPENDIX B 

We present here some more details regarding the inte- 
gration of eq. (26). Considering a spherical overdensity 
with radius R the functional form of eq.(3) becomes: 



f(R) = exp 



y 1 ^ l + w(u) 



du 



(47) 



Taking now into account that R(y) = Rty = C~ 1//3 a t t/, 
(described well in section 3.1) and differentiating the 
above integral we have the following useful formula: 

= = -s y -\i +w (R(y))]f(R(y)) -(48) 

Using now the integral of eq. (26) 

/ y[l + 3w(R(y))]f(R(y))dy = -2 / yf(R(y))dy 
+3f[l + w(R(y))]yf(R(y))dy 

and taking into account eq. (48) we can easily solve the 
integral 



/ y[l + 3w(R(y))]f(R(y))dy = -2 / yf(R(y))dy 

„2 d/(fl(y)) ,, 2 t ,r,, ( 0U > 

dy 
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